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CN ■ We calculate the effects of next-to-leading order perturbative QCD 

\ as well as of the quark transverse motion and off-shellness on the 

' Drell-Yan process cross section. By studying the s — > oo behaviour 

of the cross section in these approaches, we find that the effects of 
, quark off-shellness and intrinsic-fey parametrize those of higher twists, 

f*"""- | In particular, the off-shellness of partons generates part of the K- 

factor type corrections to the leading order cross section. Higher twist 
q . contributions to the py-spectrum of the Drell-Yan pairs are found to 

be large for presently accessible energies. The evolution of quark off- 
O shellness distribution with the hard scale is also studied. 

PACS numbers: 13.85.Qk, 12.38.Cy, 12.38.Qk 

<D '. 

> ! 1 Introduction 

X ■ 

At the future GSI-FAIR facility [1], several experiments will perform precise 
measurements of high mass lepton pair production in pp and pA collisions 
in order to determine parton distributions in the proton and in nuclei. Ex- 
periments GSI-PAX |2J and GSI-ASSIA [3] are going to study the polarized 
process, while GSI-PANDA [I] will investigate unpolarized pp — > X at s 
as low as 15 GeV 2 . 

A good theoretical understanding of the Drell-Yan process at the level 
beyond the leading order (LO) of perturbative QCD (pQCD) is necessary. 
Indeed, the effects beyond the LO pQCD are expected to be high at this 
low energy [5J [6] . These effects on the double differential cross section can be 
roughly parametrized by an overall A'-factor, giving the discrepancy between 
the LO calculations and the data. However, in case of the triple differential 
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cross section (i.e. transverse momentum distribution of dileptons), the cor- 
rections to LO pQCD cannot be parametrized by a constant factor, as we 
discuss below. 

Applied to the Drell-Yan process pp — > l + l~X, the leading order approx- 
imation of collinear pQCD predicts the correct dependence of the double 
differential Drell-Yan cross section d 2 a jdM 2 dxp on the hard scale M. Here 

%F = Pz/\Pz\max (1) 

is the Feynman variable of the lepton pair and M its invariant mass. However, 
LO pQCD fails to reproduce 

1. the magnitude of this cross section, the discrepancy being usually 
parametrized by a i^-factor; 

2. the average transverse momentum p? of the dileptons; 

3. the pr-spectrum of Drell-Yan pairs, which is given by the triple differ- 
ential cross section d 3 a/dM 2 dxpdpT- 

Experimentally observed Drell-Yan lepton pairs have non-vanishing trans- 
verse momentum p?, which can be as large as several GeV. However, in 
leading order of pQCD, the cross section is proportional to 5(pr)- Indeed, in 
collinear QCD, both the transverse momentum and the light cone energy of 
quarks inside hadrons are neglected compared to the large component of the 
quark momentum, parallel to the hadron momentum. Thus, the initial state, 
i.e. the colliding quark and antiquark to be annihilated into a lepton pair, 
has no transverse momentum. Therefore, the final state has zero transverse 
momentum, too. 

The diagrams contributing to the Drell-Yan process in addition to the 
LO parton model mechanism (a quark from one hadron annihilates with an 
antiquark from another hadron into a virtual photon) are divided into two 
classes: those of higher orders of as and those of higher twists. Calculat- 
ing next-to-leading order contributions and refitting the parton distributions 
accordingly, one can reduce the discrepancy with the data on the double- 
differential Drell-Yan cross section [TJ [H]- Moreover, the higher twists are 
vanishing in the limit of infinite energy. However, these power suppressed 
contributions can be large at realistic energies. We show in the present paper 
that the contribution of higher twists is essential for a proper description of 
the data on the triple-differential Drell-Yan cross section and propose a phe- 
nomenological model suitable to calculate these effects. We note that recently 
also H. Shimizu et al [9] have pointed out that at lower energies the dilepton 
cross section provides information on the non-perturbative dynamics. 

The lepton pairs can gain non-vanishing pt due to two possible mecha- 
nisms 

• processes of next-to-leading order (NLO) in as, 
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Figure 1: 0(as) contributions to the Drell-Yan process: (a) gluon 
Bremsstrahlung, (b) vertex correction, (c) gluon Compton scattering. Vir- 
tual photons (dashed lines) split into lepton pairs, waved lines denote gluons, 
arrows denote quarks. In each diagram, time runs from left to right. 



• quark transverse motion and off-shellness. 

In the former case, the dilepton pair recoils against an additional jet in the 
final state. An example is given by the production of a gluon besides the 
virtual photon (see Fig. [T^,) 



qq 



l+l-g. 



(2) 



Alternatively, final transverse momentum of the lepton pair can be caused 
by a non- vanishing transverse momentum of the initial state, i.e. by the 
non-collinearity of quarks inside the colliding hadrons. In this case, the 
recoil transverse momentum is carried by hadron remnants, formed by the 
"spectator" partons. In addition, the quark and gluon off-shellness can have 
a large effect for some observables, as has been recently shown in [TU] and jllj . 
The parton off-shellness arises due to interaction between the partons in one 
hadron and thus constitutes a higher twist effect. 

In the present paper, we compare the Drell-Yan process cross section and 
the transverse momentum p? of Drell-Yan pairs in both aforesaid approaches. 
We calculate in collinear pQCD the mean (pfypert, which is the part of the 



lepton p^ generated by the next-to-leading order process 



This allows 



us to study the evolution of {pj)pert with s and M. However, it turns out 
that the magnitude of the experimentally measured px cannot be described 
by NLO alone. On the other hand, the experimental data are reproduced 
much better by a model, taking into account both the intrinsic transverse 
momentum and the off-shellness of quarks in the proton, as will be shown 
in section El In this model there is no need for a i^-factor. Instead, we 
describe both shape and magnitude of double as well as triple differential 
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cross sections by extracting a quark-off-shellness from the data, which allows 
a physical interpretation. We expand the cross section in this model in 
powers of 1/s at s — > oo in order to relate the phenomenological result to 
higher twist corrections in section [3j We conclude in section HJ 



2 Effects of non-collinearity and off-shellness 
of quarks in the proton 

The transverse momentum distribution of Drell-Yan dileptons da j dM 2 dxpdp\ 
in the NLO perturbative QCD is singular at px — [T2l [To] . Neither can 
the average px of Drell-Yan pairs be reproduced in NLO: Using the method 
introduced in [12] . we calculate of dileptons in NLO perturbative QCD. 
The results are presented in Fig. [2j The px width calculated in this way is 
around 0.6 GeV at s = 1600 GeV 2 and M « 7 GeV (see Fig. [2(a)] ). This 
value is about a factor of 2 smaller than the width of the px distribution 
measured by the Fermilab experiment E866 We have to conclude that 
NLO pQCD is insufficient to describe present data on px of Drell-Yan pairs. 

In this context, we note an interesting observation: while the perturbative 
(p T ) P ert of Drell-Yan pairs increases with increasing M at high center of mass 
energy y/s (see Fig. 2(a)), it decreases with M at low s (for instance, at 



s = 32 GeV 2 relevant for the future PANDA experiment, Fig. 2(b)) due to 
phase-space limitations. 

A natural approach to generate additional px is to take into account the 
primordial transverse momentum of quarks in the proton. The primordial 
quark transverse momentum is a non-perturbative effect and, from the un- 
certainty principle, averages at > 200 MeV. On the other hand, we will show 
that the higher twist effect of the intrinsic kx on the Drell-Yan cross sec- 
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tion is of the same order as the effect of non-vanishing quark off-shellness in 
proton, which is caused by the interaction of partons of one hadron in the 
initial state. Therefore, both the intrinsic kr and quark off-shellness have to 
be taken into account for the sake of consistency. 

The data on triple differential cross section also favor a model taking into 
account both non-collinearity and off-shellness of quarks. Such a model was 
proposed and tested against the E866 data [U] in earlier publications [51 [TT]. 
The method utilizes a phenomenological transverse momentum distribution 
and an off-shellness distribution of quarks in the proton. 

This approach is based on the factorization assumption and on a second 
assumption that the soft part of the cross section can be approximated by a 
product of functions of the quark k + momentum, the transverse momentum, 
and the virtuality m 2 = k + k~ — k\. The part dependent on k + and kx is 

g q (M 2 ,£,k T ) = g(k T )q^,M 2 ), (3) 

where q(£, M 2 ) denotes parton distribution of quark flavor q at momentum 
fraction £ and scale M 2 . We use a Gaussian distribution of quark intrinsic 
transverse momentum kx 

the mean squared partonic intrinsic transverse momentum being (k T ) = AD 2 . 
We note that both the factorization ansatz (jHJ) and the Gaussian parametriza- 
tion (Tj0) are model assumptions, which are open for improvement p2], even 
though they are widely used at present [TT] . 

The quark virtuality distribution cannot be calculated from first princi- 
ples. Using analogy to many body theory, we parametrize the quark virtu- 
ality distribution as a Breit-Wigner with width T 

TT 77T + 4I 

The exact off-shell kinematics as well as the off-shell and non-collinear sub- 
process cross section (at LO in as) are used. 

The cross section of the process 'hadron A' +' hadron B' — ► l + l + X in this 
approach is [UJ 



d 3 



a 



/P POO POO pi pi 

dk±± / d%2± j dnii / dm^ / g?£i / (i£2A(mi)A(m2 
J Jo Jo Jo Jo 



dM 2 dxpdp. 



2 



X 9 Q (M , £1, kiA_)g Q (M , £ 2 , k 2 ±) , (6) 



in which not only the three-dimensional (longitudinal and transverse) mo- 
tion of partons, but also the virtualities of the active quark and anti-quark 
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are explicitly taken into account by means of a phenomenological double- 
unintegrated parton density. 

In (EJ), the following off-shell partonic cross section is used 



da 



dM 2 dxpdp. 



■ K 



T 



2\ 2 



2M 4 - M 2 (m\ - 6mim 2 + m 2 ,) - (ml - m 2 2 ) 



(ml + kf ± ) (mj + k 2 ± 
x5 \M 2 -m\- ml- ^ 2 P[ P+ - ± /V , '- + 2k 1± ■ k 2± 

r | (m 2 + k 2 ± ) (m 2 + k 2 ± 



2 



x5 [ [k 1± + k 2X ) -p 2 T ) , (7) 
with 



In©, 



3M 4 N C 8 y/(h ■ k 2 f - m 2 m 2 



(h ■ h) = te 2 p r p 2 + + K + + _ k j 2±j (9) 



(^) 2 =(^) 2 = ^-M^±^/(0 -M 2 N s (10) 

in the hadron center of mass syste 

In the scaling limit (s — > oo, s/M 2 =const), the spectral functions effec- 
tively drop out due to normalization (cf. discussion later), and the hadronic 
cross section §6§ goes to 

dM 2 dx F dp 2 T =X} $ «0 B i' a *) (dM 2 dx F ) i0 8D2 6XP (~^) ' (U) 
where 

$ ? (xx, z 2 ) = q A (x 1 )f(x 2 ) + q A (x 1 )q B \x 2 ), (12) 
( d 2 a q \ 47ra 2 e 2 Xl x 2 

\dM 2 dxF ) LO 9M 4 Xl +x 2 [ Xl%2h [ } 
and the parton momentum fractions (x\, x 2 ) are defined via 

M 2 = X!X 2 s; (14) 

x F = (x 2 — xi)/(l — x 1 x 2 ). (15) 



1 The formula ([7]) is a corrected version of the erroneous formula (39) in [llj rewritten 
in the Lorentz invariant form. Numerically, the difference between the two expressions is 
less than 1% in the kinematics studied in [TT] but might be large at xf — > 1 and pt ^> M. 
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Figure 3: Prediction for the px spectrum of Drell-Yan dileptons in our model 
at y/s = 40 GeV and xf = 0.1 as compared to the data of experiment 
E605 [20]. 

Note that the X1X2 in the denominator of the xp definition in (|14p is some- 
times omitted in textbooks, where an approximate definition xf ~ 2p z j ^fs 
is used instead of (JIJ). 

In [TT], our model was compared to the data on the triple differential 
Drell-Yan cross section d?a / dM 2 dxpdpT from experiment E866 [6] at Fermi- 
lab in pp collisions at 800 GeV incident energy. Both the slope and magnitude 
of the pt distribution of the Drell-Yan pairs were described well without the 
need for a i^-factor. In particular, the experimentally measured (p T ) is re- 
produced in this model by fitting the model parameter D (the dispersion of 
the quark intrinsic transverse momentum). At s — 1600 GeV 2 , we obtained 
D = 0.5 ± 0.18 GeV. On the other hand, the detailed shape of the distribu- 
tion turned out to be sensitive to the off-shellness, giving r = 50 — 300 MeV 
(depending on the mass bin) for this particular experiment. 

The distribution of the transverse momentum of lepton pairs produced 
in the Drell-Yan process off nuclei pA — > l + l~X also can be reproduced 
within this model. For example, in Fig. [3] the calculation for the transverse 
momentum spectrum of Drell-Yan dileptons of our model is compared to 
the data of the experiment E605 [20] on p Cu collisions at y/s = 38.8 GeV, 
xf = 0.1. The cross section plotted in Fig. [3] is 
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Figure 4: Result of LO collinear QCD for a double differential Drell-Yan 
cross section (dashed line) at yfs = 20 GeV and Xp = 0.1 as compared to 
experiment E439 [TS]. Solid line is the LO result scaled up with a factor 
K = 1.6. 



where E is given by 



The model parameters D, T used in the calculations were fitted to data on 
pp — > l + l~X in [H] and no readjustment was done for the pA case. 

In the present paper, we want to make a consistent comparison with 
NLO results. However, the triple differential cross section Ed 3 a/d 3 p, plotted 
in Fig. [HJ cannot be compared directly to the NLO result, because the latter 
is singular in every fixed order of pQCD [T31 [5]. We, therefore, apply the 
described model now to the double differential cross section d 2 a jdMdxp. The 
.fT-factor, which is needed to increase the magnitude of the LO prediction 
for the double differential Drell-Yan process cross section so that it agrees 
with the data, can be decreased from 2 to 1.1 by taking into account NLO 
processes [7]. In order to determine, what part of this LO i^-factor can be 
accounted for by the model with intrinsic kr and off-shellness of quarks, we 
compare data to the triple differential cross section (jHJ) integrated over pr- 

In Fig. HI the Drell-Yan process cross section d 2 a / dMdxF predicted at 
leading order of perturbative QCD (dashed line) is compared to the data 
of the Fermilab experiment E439 [18J on pW collision at 400 GeV incident 
energy, at xp = 0.1. The LO prediction lies below the data. The solid line 
shows the LO curve scaled up with a factor K — 1.6. The i^-factor depends 
somewhat on the parametrization of parton distributions used. We use here 
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Figure 5: Drell-Yan cross section in our model at D = 500 MeV, T = 
200 MeV compared to the data of experiment E439 [IE]. K — 1. 



the parametrization |19j . If one assumes a larger contribution of sea quarks, 
the if-factor needed to describe the data is lower. 

In order to calculate the double differential Drell-Yan process cross section 
in our model, we first find the triple differential cross section from ©• The 
double differential cross section is obtained using 

d 2 a _ r^r)™* d 3 a 2 
dM 2 dx F = J dM 2 dx F dp 2 T Pt ' ^ 

Note that the maximum transverse momentum of the Drell-Yan pair {p%) m ax 
is fixed by kinematics 

(^ (s+m :- m " )2 -4^-^. 

4s 4s 

where M\ is the minimal invariant mass of the undetected remnant. 

The data are reproduced well (see Fig. E]) with K — 1. We conclude that 
in the experimentally relevant region the .fT-factor of the double differential 
Drell-Yan cross section can be explained by two alternative scenarios: either 
as an effect of higher orders of perturbative QCD as shown in [TJ [8] or as an 
effect of non-collinearity and off-shellness of quarks in our phenomenological 
approach. The experimental cross section magnitude can be reproduced in 
NLO calculations by fitting the renormalization scale or in our model by 
fitting the parameters D and V . The latter explanation of the i^-factor has 
the advantage that it can explain also the triple differential cross section. 
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Figure 6: Drell-Yan cross section in a simplified intrinsic- fc^ approach at 
D = 50 MeV (short dash), D = 250 MeV (dash), D = 500 MeV (dash-dot) 
compared to LO pQCD (thin solid line) and the data of experiment E439 [18J. 
The solid line gives the theoretical curves multiplied by Independent K- 
factors fitted to the data. Everywhere, T = 0. 



In the following, we additionally study the relative importance of quark 
off-shellness and quark intrinsic transverse motion by comparing our result to 
that of the intrinsic-A>r approach. The intrinsic-fcr approach [21] is a limiting 
case of our model at T — > 0. The factorization assumption in this case gives 

d 4 <J ,t* . ,-> . d 2 cr(k T1 , fc T2 ) ■-/ ^ r r % / s 

g(k T1 ) ® g{k T2 ) ® TV 5(p T - kri ~ k T2 ). (20) 



dM 2 dxFdpT dM 2 dxj 

The formula is often simplified by neglecting the dependence of a on fcyi 
and kx2i for example in [13] and in PYTHIA [22J. In this case, the px spec- 
trum of Drell-Yan pairs d 3 a / 1 dM 2 dxpdp\ is also simply a Gaussian in p\. The 
cross section ( 1201) has to be integrated over the azimuthal angle of the lepton 
pair and over p\ according to (I18p . Because of the finite integration interval 
in (1181) . we do not recover the normalization of the fcr-distribution (j3J), but 
obtain a suppression that increases with D. 

The double differential Drell-Yan process cross section in the intrinsic kr 
approach with collinear sub-process cross section at three values of D is com- 
pared to the LO of pQCD and the data of the experiment E439 in Fig. [61 The 
magnitude of the measured cross section cannot be reproduced in this leading 
order intrinsic k^ approach, in which the dependence of the partonic do on 
k\T and k 2 r are neglected. Just as in LO of pQCD, scaling with an overall 
if-factor ranging from 1.6 to 4 is necessary to describe the data. It is appar- 
ent from Fig. [6] that the .fT-factor extracted from the data is D-dependent. 
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Figure 7: Drell-Yan cross section in three approaches: simplified intrinsic- 
kr with collinear sub-process cross section (dash), full intrinsic- kx (dash- 
dot), our model at Y — 225 MeV (solid). Everywhere, D = 550 MeV, 
v/i = 40 GeV, M « 7.5 GeV. 



Therefore, this scaling factor should be understood as a phenomenological 
parameter and not as a measure of higher order corrections. 

Summarizing, there are three phenomenological approaches to the Drell- 
Yan process beyond LO pQCD: 

1. our model, accounting for both intrinsic transverse momentum (D ^ 0) 
and off-shellness (r ^ 0) of quarks; 

2. intrinsic-fcT approach (D ^ 0), which is the limiting case of our model 
at T = 0; 

3. simplified intrinsic-£>r approach (r = 0), in which the primordial trans- 
verse momentum is not zero (D ^ 0), but the non-collinearity of the 
qq — > l + l~ sub-process cross section da, i.e. its dependence on ki and 
/c2, is neglected. 

We compare the effects of primordial kx, non-collinearity of da and quark off- 
shellness by plotting the triple differential Drell-Yan cross section calculated 
in the three aforesaid phenomenological approaches in Fig. [71 The simplified 
intrinsic-£/r approach gives a Gaussian for the pr-distribution (dash line). As 
we will show in the next section, the approximation of T = and collinear da 
is equivalent to restricting oneself to the leading order in the twist expansion, 
that is, in the case of the unpolarized Drell-Yan process, the expansion in 
powers of 1/M. In Fig. the importance of higher twist corrections in the 
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Drell-Yan process is illustrated by the difference between the solid and dash 
lines. 

The part of higher-twist effects incorporated in the full intrinsic- &t ap- 
proach changes the distribution considerably (cf. the dash and dash-dot 
curves in Fig. [7j). On the other hand, additional higher twist effects, mod- 
elled by quark off-shellness and given by the difference between the dash-dot 
and solid curves, are of the same order. We conclude that higher twists in 
the Drell-Yan process can be large and that we have to take into account 
both non-collinearity and off-shellness of quarks in order to model them. 

3 Twist nature of the phenomenological cor- 
rections 

In the previous section we have shown that the double differential Drell- 
Yan cross section is reproduced by a model accounting for intrinsic k? and 
off-shellness of quarks without a need for a i^-factor. In addition, the pt 
distribution of the Drell-Yan pairs can be explained in our model [UJ, but 
not in NLO of pQCD [T3l[T5] . Therefore, the effects of quark off-shellness and 
intrinsic kx do not arise solely from the diagrams of NLO pQCD. Instead, we 
will show that they parametrize higher twist processes. Some of the diagrams 
that contribute to the Drell-Yan cross section at higher twist are shown in 
Fig. [HJ Gluon radiation in the initial state and gluon exchange between the 
active parton and spectators generate intrinsic kx and virtuality of quarks in 
the proton in the Drell-Yan process. Some of these processes (for example, 
the gluon exchanges that connect factorized regions - the sub-process and 
a soft matrix element) are suppressed by powers of s in the scaling limit. 
However, the power-suppressed corrections give a sizable contribution to the 
transverse momentum spectrum of Drell-Yan pairs at finite s accessible in 
modern experiments. 

In this section, we investigate the relationship of NLO and higher twist 
corrections to those calculated in our phenomenological approach by com- 
paring their behaviour in the Drell-Yan scaling limit, in which s — > oo and 
M 2 — > oo with r = s/M 2 finite. The NLO corrections are proportional to 
as] therefore we expect only a logarithmic dependence of these corrections 
on the hard scale s ~ M 2 . On the other hand, higher twist contributions are 
suppressed in powers of s in the scaling limit. In order to determine whether 
the effects of quark virtuality and intrinsic k? are leading twist, we study 
the behaviour of the Drell-Yan cross section calculated in our model in the 
scaling limit. 

For comparison, note that the lowest order triple differential cross section 
is given by 
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From (lllj) one observes that the corrections to the pr distribution of the Drell- 
Yan pairs due to quark non-collinearity in the proton are not suppressed in 
the limit s — > oo. The model taking into account the intrinsic k? of partons 
therefore parametrizes some of the leading twist effects. This part of the 
model effects, i.e. the Gaussian smearing of p?, is related to contributions 
of the higher order diagrams of the perturbative QCD series as is shown for 
deep inelastic scattering at small Bjorken x in [Bj. 

On the other hand, the corrections to the LO cross section generated 
by the dependence of the sub-process cross section da are suppressed by 
powers of the hard scale s. Therefore, they represent part of the higher twist 
effects. To study this in more detail, we expand the cross section ([6]) in a 
series in 1/s around s = oo, keeping this time not only the leading term, as 
it has been done in (ITT)) , but all the terms that are suppressed by less than 
s 3//2 . We analyze the cross section at the specific value of x p = to make 
the formulas less bulky. 

We start from the general formula (E]). First, we expand the integrand 
of (jSJ) in 1/s. For this purpose not only da of (JZj) has to be evaluated at 
s — > oo, but also the combination of parton distributions ((3D that enters (JSJ) 
has to be Taylor expanded around (£i = \/t, £2 = V^)- The arguments of 
parton distributions £i and £2 are fixed after integrating out the ^-functions 
in As a result, the probed parton light cone momentum fractions depend 
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on quark intrinsic transverse momentum and off-shellness. After integrating 
© over £x, £ 2 and angles, the quantity $5(^1,^2) enters the hadronic cross 
section formula. Here, £1 and £2 are 



£1 



r l + 



p\j2 — ml — k 



\/rs 



s - ^( 1+ d^za +0 (i)) 



(22) 
(23) 



Keeping the first two orders in the Taylor expansion of £2) an d in the 

l/(rs)-expansion of da, we obtain: 



dM 2 dxpdp 2 -, 



o?e 2 q {l-r) 



dk\ 



=0 SvrD^r 2 ^ 3 J 
A(m 1 )A(m 2 ) exp( -^f±^ 



r?? ( n 



dm; 



(mi)r, 



PP _ 1/ 2 _ p _ p^)2 
^1^2 "'l ^ 



G\{t, 



(Pt 2 1 2 

— I y - m 2 - k 2 



+ T q (r 



2 7 2 
m x — k-t 



TS „, 

I h r (mi, m 2 

8 16 V ' ' 



(24) 



where 

mm 
max 

(mi) 
(w 2 ) max: 
F(m 1 ,m 2 ) 



= (p T - k 2 ) 2 ; 
= ( PT + k 2 ) 2 ; 



= ^rs+ P 2 T /2-kl; 
= \ ts+p 2 t /2-% 



(25) 
(26) 

(27) 
(28) 



2{m\ + m^(r/8 + g/6 - + (mi + m 2 ) 2 r/6 



-m 



?r/6 + (m? - ^)liV^/6 - m?rV^/(6|i) 



+ T -{m\ -m\ + k\- k 2 2 )^j 

t%t) = ^(v^v^); 



and 



G\{t) 



GI{t) 



d$ q {xi,x 2 ) 



dxi 
9$ 9 (xi,x 2 ) 



dx-y 



{x 1 =^t,x 2 =^/t) 



{x\=sJ~T ,X2 = \/t) 



(29) 
(30) 

(31) 
(32) 



are the derivatives of the parton distribution product around (y/r, \^t). Note 
that the term in square brackets in eq. ( 1241 is not symmetric in k\ and k 2 . 



14 



This is due to the different integration boundaries present in the k\ and k 2 
integrations. Several (^-functions are evaluated to get from (jSj) to ([21]). In 
the course of these calculations, it turned out to be useful to treat k\± and 
&2J_ differently. 

To further investigate the dependence of the integral (12111 on s, we have 
to specify the quark spectral function. Indeed, the integration variables m% 
and m 2 at s — > oo can be arbitrarily big, as can be seen from (1271) and (I28p . 
Therefore, only after the integration over rrii has been performed can we 
judge whether any off-shellness generated term is sub-leading in s and how 
much it is suppressed. On the other hand, the integration over rrii provides 
additional terms ~ k 2 /s, interconnecting the off-shellness and intrinsic 
effects. 

In the following, we perform the analytical integration of (j21J), assuming 
different functional forms for the spectral function A(m): 

1. a Dirac delta-function 5 (to), 

2. a Breit-Wigner function (Lorentz distribution) with a constant param- 
eter T, see (j^l). 

In the former case, the model reduces to the intrinsic-A/p approach. Inte- 
grations over rrii drop out, while the remaining integrals over k\ and k\ can 
be done analytically via Bessel functions. As the result, one finds the leading 
term ( TTTT) plus \j(rs) suppressed contributions. 

Let us now consider the second, more general, case. The cross section for 
A{m) = 5(m) is the limiting case of the formulas given below for a Breit- 
Wigner distribution (JHJ) at T = 0. Inserting the spectral function (jSJ) into 
(I24p . performing all the integrations and keeping only the first few leading 
terms in 1/M, we obtain (note that M 2 = ts — > oo, as s — > oo) 



d 3 a q 



dM^dxpd'pTp 



_r^ exp (-ft?(^), [T9(r) 

+ {^V)-^(G!(r) + Gl(r))}^ 

, _ T , [T) [ ^i + jD 2M i 



3 W U J J M 2 



(33) 



At leading twist, the Gaussian distribution of pt (111!) is recovered. How- 
ever, it is modified by the higher twists, suppressed in the limit s — > oo, but 
substantial at finite s accessible in experiment. The term proportional to 
1/M = 1/y/rs is pr-independent and leads to an overall enhancement of the 
cross section, while the pr-dependent terms proportional to 1/M 2 addition- 
ally modify the shape of the pr distribution. 
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The contribution of the off-shellness of quarks to (1331) is given by the 
summands proportional to T. It is suppressed by powers of M and vanishes 
in the intrinsic hr approach, in which T = 0. Thus, the model, which ad- 
ditionally accounts for quark off-shellness, parametrizes more higher twist 
effects than the intrinsic-/cT approach alone. 

It is interesting that the effects due to the finite quark width T appear 
in the expansion at odd powers of 1/M in contrast to those due to the 
intrinsic- The first T-dependent correction is proportional to 1/M = 
1 / y/rs. Therefore, the corrections due to the virtuality of quarks seem to 
have a non-analytical dependence on s as (rs) -1 / 2 . In order to preserve 
analyticity of the cross section we have to assume that the quark spectral 
function width T has a particular scaling behavior at large hard scale of the 
probe M = y/rs: 

r(M) ~ i as M -> oo. (34) 

Then, in ( 1331) . the terms proportional to Y/M and the terms proportional to 
1/M 2 together constitute the dominant higher twist correction to the leading 
result (jlip in the scaling limit. 

We expect the formula fl3Hj) to give a good approximation to the Drell- 
Yan cross section ([6]) at large finite M and s. In order to illustrate this, we 
compare the result of the exact calculations, i.e. the numerical integration of 
(jEj), to the leading twist approximation ( JTT1) and to the next-to- leading twist 
result (l3"3"j) in two regimes: 

• at M « 7 GeV and s = 1600 GeV 2 , see Fig. [HJ 

• at M = 1 GeV and s = 30.25 GeV 2 relevant for FAIR [T], see Fig. QH 

As expected, the sum of leading and next-to-leading terms of the power 
series ( 13H]) reproduces the full calculations quite well at M as high as 7 GeV. 
The approximate cross section has the same average magnitude and slope. 
Therefore, it is dominating the i^-factor type corrections to the leading twist 
cross section. Only the bend of the cross section at low pt, which is seen 
in the full calculations and in the data (Fig. [3]), is not reproduced at the 
next-to-leading twist. 

From the Fig. [TTJ], one sees that our model predicts the higher twist effects 
to be very large at low M and s. The discrepancy between approximate and 
exact Drell-Yan cross sections is large in this regime, too, especially at low 
Pt- We conclude that one needs to go beyond the next-to-leading twist at 
this low M and s. In this region, our model becomes indispensable, because 
it effectively sums higher orders and higher twists. 

4 Conclusions 

We have analyzed the double differential Drell-Yan cross section d 2 a/dM 2 dxF 
and the pt distribution of the Drell-Yan dileptons d 3 a/dM 2 dxpdp T in two 
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P T , GeV 

Figure 9: Cross section of pp — > l + l~ X at M ~ 7 GeV, s = 1600 GeV 2 in our 
model (solid), in the leading order in 1/M 2 (dash), up to the next-to-leading 
order in 1/M 2 expansion (dash-dot). D = 650 MeV, T = 225 MeV, x F = 0. 
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0.0 0.5 1.0 1.5 2.0 2.5 

P T , GeV 

Figure 10: Cross section of pp -> l + l~X at M = 1 GeV and s = 30.25 GeV 2 
in our model (solid), in the leading order in 1/M 2 (dash), up to the next-to- 
leading order in 1/M 2 expansion (dash-dot). D = 600 MeV, T = 250 MeV, 
xf = 0. 
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alternative approaches: collinear perturbative QCD at next-to-leading order 
and a model, which makes use of phenomenological distributions for and 
off-shellness of quarks in the proton. 

We find that the transverse momentum spectrum of the Drell-Yan pairs 
at the next-to-leading order pQCD disagrees with experiment both quanti- 
tatively and qualitatively. In contrast, we find that the phenomenological 
model with off-shell non-collinear partons successfully describes both the 
double differential Drell-Yan cross section and the pt spectrum of Drell-Yan 
pairs without the need of a K factor. 

The analysis of the Drell-Yan process cross section in our model in the 
Drell-Yan scaling limit has shown that the phenomenological model parametrizes 
higher twist effects. Higher twist contributions were up to date usually con- 
sidered to be small, because they are suppressed by powers of the hard scale. 
As a rule, they are neglected in pQCD calculations. However, the power 
suppressed effect can be large at realistic energies. 

We have found that the intrinsic transverse momentum of quarks gener- 
ates both leading twist and 1/(ts) = 1/M 2 suppressed effects. This is in line 
with our analysis of section [H which has shown that only part of observed 
{pj) can be explained by NLO effects. In addition, we have shown that next- 
to-leading twist corrections due to quark off-shellness lead to an overall cross 
section suppression and are therefore responsible for a part of the K-factor 
type discrepancy between the leading order pQCD and data. Although we 
have shown this only for a Breit-Wigner parametrization of the spectral func- 
tion, we do not expect these results to depend on its special form, since all 
calculated cross sections involve only integrals over the spectral functions. 

If a Breit-Wigner parametrization for a quark spectra function is used, 
the next-to-leading contribution is proportional to Y/y/rs. This leads us to 
suggest that the quark spectral function width Y scales as Y(M) ~ 1/M at 
large hard scale M = \frs. 

The formula that we obtained for the Drell-Yan cross section at the next- 
to-leading twist level can be very useful for applications, for example, in an 
event generator. Indeed, it requires no numerical integration, while providing 
a good approximation to the full calculations at M > 5 GeV. However, at 
M < 5 GeV, one has to go beyond the next-to-leading order in the power 
series and use the formulas of [TT] . 

The results show that the higher twist corrections to high energy processes 
can be large. Therefore, a detailed study and modelling of these effects is 
necessary, if one hopes to reliably extract quark and gluon properties from 
hadron scattering data. 
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